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INTRODUCTION 


This  report  constitutes  a  Semi-Annual  report  for 
contract  N00014-73-C-0270  administered  by  the  office  of 


Naval  Research.  The  research  is  being  carried  out  at 

3 

Systems,  Science  and  Software.  (S  ) 


The  research  performed  under  this  contract  is  a 
continuation  of  work  begun  under  a  previous  contract  at 
S  to  investigate  methods  for  determining  the  sensitivity 
of  the  results  of  complex  calculations  to  uncertainties 
in  the  modeling  parameters  entering  the  calculations. 

The  emphasis  to  date  has  been  concentrated  on  non- 
equilibrium  chemical  kinetics  systems,  to  determine  the 
sensitivity  of  calculated  concentrations  to  experimental 
uncertainties  in  the  rate  coefficients. 


In  the  previous  contract,  a  method  called  the 
Fourier  Amplitude  Sensitivity  Test  (FAST)  was  developed 
and  applied  to  two  simple  chemical  systems.  It  was  shown 
that  the  method  does  provide  a  measure  of  the  relative 
•contribution  from  each  rate  constant  uncertainty  to  the 
uncertainty  of  the  calculated  concentrations. 

The  present  research  effort  is  just  getting  under 
way,  and  consists  in  applying  FAST  to  two  particular 
chemical  kinetics  problems  of  interest  to  the  high 
altitude  fireball  community.  Results  from  this  work 
will  be  available  soon. 
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The  work  that  was  done  in  the  previous  contract  has 
recently  been  written  up  in  two  companion  papers  which 
have  been  accepted  by  the  Journal  of  Chemical  Physics  for 
publication.  Preprints  of  these  two  papers  are  included 
as  the  remainder  of  this  report. 


'"'  '  '  ‘  ’"  '''  "  "’"V'1'  '"■-  '  ’’  "  '  '  "■“  "  ”  '  '  •'■*,!•«■_>•■■  •■■  *•  "'!  '••'•  ’’-  " 

acr^rs»wK-?^iii)««ww»ii.«f ,y',»  mvw'g-aB-’wiwwKwm^^ 


A  STUDY  OF  THE  SENSITIVITY  OF  COUPLED  REACTION 
SYSTEMS  TO  UNCERTAINTIES  IN  RATE  COEFFICIENTS 
I.  THEORY  * 


R.  I.  Cukiert  C,  M.  Fortuin  and  K.  E.  Shuler 
Department  of  Chemistry 
University  of  California,  San  Diego 
La  Jolla,  Cal .  92037 


and 


A.  G.  Petschek**  and  J.  H.  Schrnbly 
Systems,  Science  and  Software 
P.  0.  Box  1620 
La  Jolla,  Cal .  92037 


*  Supported  in  part  by  the  Advanced  Projects  Agency  of  the  Department  of 
Defense,  monitored  by  the  U.S.  Office  of  Naval  Research  under  Contracts 
No.  N00014-69-A-0200-601 8  and  N00014-71 -C-J347( POOUD2) . 


I 


I 


S 


1 


ABSTRACT 


A  method  has  been  developed  to  investigate  the  sensitivity  of  the  solutions 
of  large  sets  of  coupled,  non-linear  rate  equations  to  uncertainties  in  the  rate 
coefficients.  This  method  is  based  on  varying  all  the  rate  coefficients  simul¬ 
taneously  through  the  introduction  of  a  parameter  in  such  a  way  that  the  output 
concentrations  become  periodic  functions  of  this  parameter  at  any  given  time  t. 
The  concentration  of  the  chemical  species  are  then  Fourier  analysed  at  time  t. 

We  show  via  an  application  of  Weyl's  ergodic  theorem  that  a  subset  of  the  Fourier 
coefficients  is  related  lo  (^r~y>  The  rate  of  change  of  the  concentration  of 
species  i  with  respect  to  the  rate  constant  for  reaction  %  averaged  over  the 
uncertainties  of  all  the  other  rate  coefficients.  Thus  a  large  Fourier  coeffi¬ 
cient  corresponds  to  a  large  sensitivity,  a  small  Fourier  coefficient  corresponds 
to  a  small  sensitivity.  The  amount  of  numerical  integration  required  to  calcu¬ 
late  these  Fourier  coefficients  is  considerably  less  than  that  required  in  tests 
of  sensitivity  where  one  varies  one  rate  coefficient  at  a  time,  while  holding 
all  others  fixed.  The  Fourier  method  developed  in  this  paper  is  not  limited 
to  chemical  rate  equation,  but  can  be  applied  to  the  study  of  the  sensitivity 
of  any  large  system  of  coupled,  non-linear  differential  equations  with  respect 
to  the  uncertainties  in  the  modeling  parameters. 
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I.  Introduction 


Sets  of  coupled,  non-linear  rate  equations  arise  in  a  number  of  disciplines. 
A  classic  example  is  that  of  chemical  rate  equations.  In  ne  study  of  combustion, 
air-pollution,  upper  atmosphere  phenomena  and  chemical  lasers  as  many  as  100 
coupled  rate  equations  involving  some  50  separate  species  may  be  needed  to 
account  for  the  properties  of  such  systems.  One  is  then  faced  with  the  problem 
of  solving  a  large  set  of  coupled,  non-linear  differential  equations  of  the  form 


=  Ffcjik} 


(l.D 


when  c  is  a  vector  of  concentrations,  { k)  a  set  of  rate  coefficients  and  ?  some 
complicated  function  of  the  concentrations.  While  one  may  prove  existence  and 
stability  theorems  for  the  equilibrium  point(1),  the  only  way  to  solve  these 
equations,  i.e.  to  obtain  all  the  species  concentrations  as  functions  of  time, 
is  through  the  use  of  a  high-speed  computer. 


Unfortunately,  the  rate  coefficients  (or  cross  sections)  for  many  reactions 
are  not  known  with  high  accuracy  and  indeed  may  be  uncertain  by  one  or  more 
orders  of  magnitude.  This  gives  rise  to  the  very  important  problem  of 
sensitivity"  which  may  be  defined  as  the  effect  of  uncertainties  in  the  rate 
coefficients  on  the  calculated  concentrations  of  all  the  various  intermediate 
and  product  species.  The  uncertainty  in  the  rate  coefficients  of  certain 
important  reactions  in  the  reaction  scheme  may  have  a  significant  effect  on 
the  .output  function  (for  instance,  concentration  at  time  t),  while  uncertainties 
of  the  same  magnitude  in  rate  coefficients  of  "unimportant"  reactions  in  the 
reaction  scheme  may  hardly  effect  the  output  function.  The  reliability  of  the 
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output  numbers  clearly  cannot  be  established  without  knowledge  of  the  sensitivity 


of  the  output  data  to  the  uncertainty  in  the  input  parameters. 


The  problem  is  to  find  a  practical  method  of  determining  the  effects  of 


the  uncertainties  in  the  rate  coefficients  on  the  solutions  of  the  rate  equations 


Since  we  are  interested  in  the  situation  where  the  uncertainties  in  the  rate 


coefficient  may  be  orders  of  magnitude,  linearization  schemes  are  not  appropriate, 


A  "brute  force'  method  of  investigating  the  sensitivity  is  not  feasible  as  can 


readily  be  seen  from  the  following  example.  Suppose  we  have  a  reaction  scheme 
which  has  n  coupled  reactions  involving  m  different  chemical  species.  Let  us 


furthermore  assume  that  we  wish  to  calculate  the  concentrations  of  the  m  species 


at  some  time  t  for  z  different  values  each  of  the  2n  rate  coefficients.  If  we 


now  change  one  rate  coefficient  at  a  time  while  holding  all  the  others  fixed, 
we  would  have  to  carry  out  z2n  integrations  of  the  rate  equation  (1.1)  to  time  t. 
For  the  m  different  species,  this  procedure  will  give  rise  to  a  print-out  of 
m(z)  concentrations.  If  we  know  to  a  good  accuracy  the  equilibrium  constants 
for  all  the  reactions  and  apply  the  principle  of  detailed  balance,  the  number 
of  independent  rate  coefficients  will  be  reduced  to  n  and  we  would  have  to  carry 
out  zn  integrations  up  to  time  t  for  each  species  m.  In  either  case,  for  n 
large,  it  is  obvious  that  the  time  and  expense  involved  in  such  an  analysis  of 
sensitivity  is  prohibitive  and  the  print-outs  so  numerous  as  to  defy  a  simple 
analysis  of  the  results.  Clearly,  one  needs  to  devise  some  more  powerful 


method  for  the  study  of  sensitivity. 


Our  approach  to  this  problem  is  to  ask  for  a  less  detailed  description  of 
the  effect  of  rate  coefficient  uncertainty  on  the  output  function  at  any  giv  »n 


time.  We  vary  all  the  rate  coefficients  simultaneously  so  as  to  explore  the 


entire  space  of  uncertainties  in  the  rate  coefficient  set  { k} .  As  will  be  seen 
below,  this  turns  out  to  be  equivalent  to  varying  a  single  rate  coefficient  and 
then  averaging  tne  attendant  concentration  changes  over  the  uncertainties  of  all 
the  other  rate  coefficients,  where  we  express  this  uncertainty  in  terms  of  a 
probability  distribution.  Our  approach  is  thus  related  to  a  "mean  field"  theory 
where  we  represent  the  fluctuations  of  the  field  by  the  uncertainties  in  the 
(n-1)  rate  coefficients  over  which  we  average. 

To  carry  out  this  program,  we  relate  each  rate  coefficient  k.  to  a  frequency 
introduce  a  parameter  s  which  simultaneously  varies  all  the  rate  coef¬ 
ficients  in  such  a  way  that  the  concentrations  at  a  fixed  time  become  periodic 
functions  of  s.  The  concentrations  can  then  be  Fourier  analyzed.  We  then  show 
that  a  certain  subset  of  these  Fourier  coefficients  can  be  related  to  the  first 
partial  derivative  of  tne  concentration  c.  of  species  i  with  respect  to  a  rate 
coefficient  k£  averaged  over  the  uncertainties  of  all  the  other  rate  coefficients. 
A  large  value  of  the  Fourier  coefficient  then  shows  that  /~Ais  large, 

a  Vk2/ 

i.e.  the  effect  of  a  change  in  the  a'th  rate  coefficient  on  the  concentration 


of  species  i  is  significant.  Conversely,  a  small  Fourier  coefficient  A^ 

/Sci\  .  “j 

indicates  tnat/^-y  is  small,  i.e".  the  effect  of  the  variations  of  the  j'th 

rate  coefficient  on  the  concentration  of  species  i  is  small.  In  order  to  calcu¬ 
late  these  Fourier  coefficients,  we  must  integrate  the  rate  equations  numerically 
up  to  the  desired  time  for  each  value  of  the  parameter  s.  The  number  of  s  values 
which  we  include  in  our  parameter  set  determines  the  accuracy  to  which  we  can 


calculate  the  Fourier  coefficient;  the  larger  the  set  of  s  values,  the  more 
accurate  the  determination  of  the  Fourier  coefficients. 

Since  we  must  still  perform  numerical  integrations  of  the  rate  equations 
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0-1 )»  the  question  arises  why  this  method  is  to  be  preferred  to  the  more 
direct  method  of  varying  each  rate  coefficient  separately  while  keeping  all 
others  fixed.  As  will  be  shown  in  paper  II,  which  deals  primarily  with  the 
computer  calculations,  the  number  of  integrations  required  by  the  Fourier 
method  is  0(n  )  where  r  is  a  small  integer  (r<10)  which  depends  upon  the  choice 
of  the  frequencies  i  =  l,2,...,n.  It  can  readily  be  verified  tiiat  for  ri, 
z  large,  n  <(z)  .  The  computational  economy  of  this  method  of  analysis  thus 
becomes  more  pronounced  as  n,  the  number  of  reactions,  increases.  The  reason 

for  this  reduction  in  the  number  of  required  integrations  up  to  time  t  is  to 
be  found  in  the  fact  that  in  the  Fourier  method  we  sample  the  c(k)  space  at 
a  set  of  points  determined  by  the  values  of  the  set  {s}  and  the  vector  £, 
whereas  the  "brute  force"  method  involves  neithe>  sampling  nor  the  simultaneous 
variation  of  the  set  of  rate  coefficients  { k}  and  thus  requires  many 

more  integrations  of  the  r^te  equations.  As  will  be  clear  from  the  body  of 
the  paper,  this  sampling  in  a  certain  sense  corresponds  to  the  averaging  over 
all  the  rate  coefficients.  The  reduction  in  the  number  of  required  integrations 
is  thus  intimately  related  to  the  use  of  a  "mean  field"  theory. 

Our  results  are  presented  in  two  papers,  I  and  II.  In  paper  I,  we  present 
the  theoretical  basis  of  our  method  without  explicit  reference  to  the  verifying 
computer  experiments.  In  paper  II,  we  present  the  results  of  our  computer 

calculations  as  well  as  a  detailed  discussion  of  the  problems  involved  in  such 
calculations. 

It  should  be  pointed  out  that  the  utility  of  this  Fourier  analysis  method 
of  testing  sensitivity  extends  beyond  the  confines  of  chemical  kinetics  and 
beyond  the  confines  of  differential  equations.  Large  sets  of  coupled,  non¬ 
linear  equations  are  used  in  many  fields  such  as  economics,  population  dynamics, 
weather  forecasting,  systems  analysis,  operations  research,  etc.  for  modeling 
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and  predictive  purposes.  As  has  been  pointed  out  by  a  number  of  investigators, 
it  is  important  that  sensitivity  tests  be  carried  out  on  such  systems  to 
identify  the  critical  parameters  and  to  validate  the  applicability  of  the 
models.  The  method  developed  here  can  be  applied  to  any  set  of  equations, 
differential,  integral,  algebraic,  etc.,  which  yield  an  output  as  a  complex 
function  of  many  parameters.  In  fact,  we  have  used  an  analytic  function  of 
several  variables  to  test  some  of  our  ideas.  We  plan  to  apply  this  Fourier 
method  to  the  analysis  of  other  complex  systems  in  the  near  future. 
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II .  Fourier  Analysi s 


The  rate  equation  for  coupled  chemical  reactions  which  obey  the  law  cl 


mass  action  can  be  written  in  the  form 


(2.1) 


( ^  ^ )  •  •  •  ,m) 

where  c  ((k, ...  .Jcn;t)  is  the  concentration  of  species  i  at  time  t,  or  .  .  ..  v' 

'S  the  stoichiometric  coefficient  of  the  species  1  in  reaction r,  with  r=l . n 

labeling  the  different  reactions  in  the  reaction  system,  and  where  k  (k  )  is 
the  forward  (backward)  rate  coefficient  for  reaction  r.  The  coefficients  y  f 
and  vr>,.are  non-negative  integers  so  that  the  stoichiometric  coefficient 
vrj  defined  above  can  be  positive,  negative  or  zero.  From  the  form  of  the  rate 
equations  (2.1)  one  can  prove  that  the  concentrations  are  bounded  and  non-negative, 

that  they  are  continuous  functions  of  the  time  and  the  rate  coefficients  and 
that  an  equilibrium  nninr  pvktc  0) 


that  an  equilibrium  point  exists.^ 

In  order  to  determine  the  effects  of  uncertainties  in  the  rate  coefficients 

kir  on  the  concentrations  Cj  at  time  t  a  systematic  method  for  varying  the  k's 
must  be  developed.  We  define 


kf  =  k!°U 


i  i 


0=1 . rr>) 


(2.2) 


Ui"fi(Sln“i-S)  (2.3) 

where  k!o)is  the  "best  value"  of  the  rate  coefficient  (i.e.  the  one  which  the 
investigator  believes  to  be  the  best  available  value. based  on  experiments  or 
calculations),  the  "frequency"  „.  is  a  positive  integer,  s  is  a  parameter  and  f 
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is  a  function  to  be  determined.  The  introduction  of  the  parameter  u.  and  the 
form  of  Eq.  (2.2)  permits  one  to  effect  readily  oroer  of  magnitude  changes  in 
the  rate  coefficients.  The  form  of  equation  (2.3)  permits  one  to  vary  simul¬ 
taneously  all  the  rate  coefficients  by  varying  tne  parameter  s.  The  rapidity 
of  the  variation  is  determined  for  each  k.  by  the  magnitude  of  the  w.'s.  The 

ui  are  Chosen  t0  be  P°sitive  integers  in  order  that  the  concentrations  at  a 
fixed  time  t  become  periodic  functions  of  s  with  period  2tt: 

ci(s;t)  =  c.(s  +  2-n ;  t)  (2 

In  the  development  to  follow  we  suppress  this  dependence  of  the  concentrations 
on  time,  it  being  understood  that  the  analysis  is  carried  out  for  a  fixed  time  t. 

The  concentration  as  a  function  of  s  describes  a  closed  path.  That  is, 
for  each  1  and  for  every  value  of  s  we  obtain  a  point  in  t  space  with  value 
c.jOO;  as  s  changes  by  2*  we  return  to  the  same  values  of  t  and  u,  [see  eqs. 

(2.2)  and  (2.3)],  and,  from  equation  (2,4),  to  the  same  value  c.(s).  Since 
ci(s)  is  pe»  iodic  on  2tt  we  may  expand  it  in  a  Fourier  series 

/\(  i)  “ 

Cl'(s)  =  ~2  +  “  (Ar^  S1'n  rs  +  cos  rs)  •  (2.5) 

In  the  analysis  presented  below  we  are  interested  only  in  the  Fourier  sine 

coefficients  that  correspond  to  the  original  input  frequencies  a,.,  i.e.  the 
coefficients  given  by 


A<’'>  .  i 
“l  * 


sin  co^s  ds 


(2.6) 


£  =  1 ,2,. . . ,n 

We  now  wish  to  relate  the  Fourier  coefficients  to  the  effect  of  the 
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variation  of  the  rate  coefficients  k^  and  the  concentrations  c . .  In  order  to 


do  this  it  is  iiecessary  to  relate  the  s-  space  integral  of  Eq.  (2.6)  to  a  n- 


dimensional  integral  over  the  entire  u  space.  It  can  be  shown  that  if  the 


entire  u  space  is  covered  densely  as  s  is  varied,  then  the  integrals  over  the 


s-  space  and  the  u  space  yield  identical  resalts.  This  is  just  a  statement 


of  the  ergodic  theorem  which  permits  one  to  equate  time  and  phase  space  averages 


in  statistical  mechanics.  If  the  frequencies  u.  in  Eq.  (2.3)  were  chosen  to 


be  incommensurate,  then  the  function  c.(s)  i.e.  the  concentration  as  a  function 


of  the  parameter  s,  would  be  an  almost  periodic  function.  This  implies  that 
the  path  in  the  space  Rn(3),  where  the  0's  are  defined  by 


0^  =  u  •  s  (mod  2tt) 


(2.7) 


returns  arbitrarily  close  to  any  initial  point  as  s— One  could  then  equate 
the  s-  space  integral  with  the  integral  over  the  n-  dimensional  9  space  as  was 
first  proved  by  V/eyl . ^  ^  However,  the  use  of  an  incommensurate  set  of  frequencies 
“i  in  Ecl-  (2-3)  would  require  that  the  numerical  evaluation  of  the  Fourier 
coefficients  of  Eq.  (2.6)  be  carried  out  over  an  infinite  period.  This 


clearly  is  not  a  feasible  procedure  on  any  known  computer. 


It  is  for  this  reason  that  we  introduce  integer  frequencies  u.  in  Eq.  (2.3) 
Their  use  leads  to  a  finite  period  (0  to  2tj)  analysis  which  can  readily  be 


handled  on  a  computer.  Tnis  will  give  rise  to  an  error  in  the  analysis  since 


now  the  "phase  point"  will  no  longer  densely  cover  the  u  space  as  s  is  varied 
and  the  s-space  and  u  space  integrals  therefore  do  not  yield  identical  results. 


Let  us  for  tne  moment  ignore  this  error  and  use  Weyl's  theorem  for  oi 
periodic  function  c(s).  We  write 


"T n  dej  ci<e, . 

»  O  ■'n  i—l 


(2.8) 


/o  j=l 


V-  Mi’tsnb  is  ■&(.***  U 


Tin's  expression  can  be  put  in  a  more  suggestive  form  as  follows.  The  £ '  th 
integral  is  rewritten  by  integration  by  parts  as 


L  f  d0£  c(9i'--.en)  sin  0£  =  — 
•no 


dS£  Hr cos  st 

X 


(2.9) 


where  the  boundary  term  has  vanished  by  the  periodicity  of  the  integrand.  The 
use  of  Eq.(2.9)then  permits  us  to  rev/rite  Eq .  (2 . 8 }  as 


l(1)  -  2 


C 271  f271  n 

a  "in 


deJ  5^  ‘“s  et 


(2.10) 


The  6  space  integral  must  now  be  related  to  the  u  space  integral.  The  transfor¬ 
mation  is  (see  Eqs.  (2.3)  and  (2.7)) 
df.(sin  0.) 


duJ  =  -  a  s'l'n  6J  cos  ej  dl)j 

j 

where  we  require  that  f.  be  a  monotonic  function  of  its  argument.  Then 

J 


(2.11) 


,  ,  fM+1)  rfn(+i ) 

p.w--2-\  ...( 

..  n  I  ; 


£  n 


^H)  fn("D 


adc. 

duj cos  \  af 

~n  df.  (sin  9~ 

j!  cos  0-  : — ~- 

j/£  J  5  51 "  9j 


(2.12) 


In  order  to  obtain  a  definite  relation  between  A'1  ^  and  dc./Su  we  must  choose 

U)£  1  £ 

some  particular  form  for  the  function  f j (sin  0 j ) .  There  are  several  choices  which 
lead  to  useful  expressions  for  the  A^;  as  will  be  seen  below,  a  particularly 


advantageous  choice  is  to  set 


aVsin0j) 


b  sin 


cos20.  =  —  , 
J  3J 


(2.13) 


where  aj  is  a  parameter.  The  use  of  Eq.  (2.13)  in  Eq .  (2.12)  leads  to 
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k(1)  -  2 


/•■+«  /•■*«  n 

n  J  •')  TT 

*,-“  ^  j=l 


£  a£n 


du .  a  .  cos  0 .  — - 
J  J  j  du£ 


(2.14) 


Integration  of  Eq.  (2.13)  yields 

f.fc-in  n  \  _  „  _  1  fl  +  Sin  6.' 

ijisin  u- )  -  u,  -  £n  = - : — rtj. 

J  J  J  2a •  1  -  sin  0. 

j  jj 

and  also  indicates  the  range  of  integration  in  Eq.  (2.14). 
as  a  function  of  Uj  from  Eq.  (2.15)  gives  for  Eq.  (2.14) 

/•-foo  --foo  n 


(2.15) 


Expressing  cos  0. 

J 


Since 


Tf)  „ 

2 

\ 

,  n  . 

V 

••-fco 

r-K°  n 

Tf 

L". 

1 1 

-X  j=l 

/•+O0  *+«  n 

;  -J  it  o 

7T  —X)  J^o  -j  r:  I  J 


cosh  a.u  du 
J  J  ^ 


du.  - i - - -  =  T 

J  cosh  a.u. 

J  J 


we  obtain  as  our  final  result 

A“)’ =  <iv3u*> 

where  the  bracket  in  Eq.  (2.18)  is  defined  by 

r°°  n 

I  •••  J  TT  P<ui ;  a  )  Y(u . uj  du. 

(Y(u  . u  )>  =  ~ _ J-“  J=1  J  J  1  n  J 

l  n  /•+”  /-+c°  n  ‘  - 

'  TT  p(u.;  a.)  du. 

•'-oo  J  —  1  J  J  0 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


* 

p(u,;a.)  =  — L - 

J  J  cosh  a.u. 

J  3 


(2.20) 


The  function  p(Uj ;a^ )  can  be  interpreted  as  a  distribution  function  in  2  space 
which  weights  the  uncertainty  in  the  rate  coefficients.  Equation  (2.18)  is  the 


desired  relation  between  the  Fourier  coefficient  ,V  1 '  and  the  chanqe  in  the 

concentration  of  species  i  witn  a  change  in  the  L’tii  rate  coefficient,  Dc./Du  , 

1  *• 

averaged  over  the  changes  of  all  the  other  rate  coefficients. 


The  particular  form  of  f.  (sin  0.)  that  we  have  chosen,  Eq.  (2.15),  has 

J  J  '  * 


lead  to  a  weight  function  p(u.;a.)  for  each  rate  coefficient,  Eq.  (2.20),  which 

•J  J 

has  several  convenient  properties.  As  a  function  of  u.,  the  function  p(u.;a.) 

J  J  J 


is  symmetrically  bell  shaped  about  u.  =  Q  corresponding  to  k.=k(°b  the  "best" 

J  0  J 


value  of  the  rate  coefficient  k..  As  a.,  which  is  a  parameter  at  our  disposal, 

J  J 


is  increased,  the  weight  function  p(u.;a.)  narrows  about  u.-O;  this  corresponds 

J  J  J 


to  a  decreased  spread  in  the  values  of  the  rate  constant  k.,  i.e.  it  corresponds 

to  a  narrower  range  of  the  uncertainty  of  k.  around  k(0/ .  In  the  limit  as  a-r- 

J  J  J 

one  obtains 


lim  p  ( u  . ;  a  . )  =  tt6(u.) 
aj— J  J  J 


(2.21) 


where  6 (x )  is  the  delta  function  of  argument  x  which  implies  that  the  rate  con¬ 


stant  is  known  with  certainty  to  be  kj0^.  When  some  information  is  available 


on  the  spread  and  distribution 
,(o) 


of  the  rate  coefficient  k.  about 

J 


its  "best"  value  k^  ,  one  can  determine  the  a.  through  the  standard  deviation 

j  j 


of  the  values  in  u  space 


<u?>  .  I 

J  7T 


/•« 

/  - 

•^-00 


3 

cosh  a.u. 

J  J 


1 


2_  ■  I 

4  *5 


(2.22) 


from  which  it  follows  that 


a,  =  i  1 


j  2  (u2)1/2 


(2.23) 
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Ttic  parameter  a.  thus  permits  us  to  introduce  explicitly  the  spread  of  uncertain- 

J 


ties  in  the  values  of  the  rate  coefficients  into  our  analysis. 


To  calculate  the  Fourier  coefficients  A^  we  must  first  choose  a  suitable 


w. 


set  of  w j ' s  and  a^.'s.  What  constitutes  a  "suitable"  set  of  id's  will  be  discussed 


in  the  next  section  in  conjunction  with  the  error  analysis  and  in  more  detail 

in  paper  II.  The  set  of  u-'s  defines  the  u.'s  as  a  function  of  s  according  to 

J  J 


Eqs.  (2.3),  (2.7)  and  (2.15).  For  each  value  of  s  one  obtains  a  value  of  the 


concentration  c  .  (k^ , . . .  ,kn) .  In  principle  one  can  then  compute  the  to  any 
desired  accuracy  by  taking  enough  values  of  s, 


-  fiiiy. 


m 


,  q  =  1 ,2, . . . ,m, 


(2.24) 


where  ni  is  some  integer. 

It  is  important  to  point  out  that  our  main  interest  is  in  identifying  those 


rate  coefficients  whose  variation  significantly  effects  the  concentration  c. 


of  a  species  i  at  time  t,  and  those  rate  coefficients  whose  variation  has,  only 


a  minimal  effect  on  the  species  concentration  c..  Thus  if  one  of  the  Fourier 


M) 


coefficients,  say  Av /,  is  one  or  more  orders  r>f  magnitude  larger  than  all 


other  coefficients  Av  i  =  l,2,...,n,  (t/j),  than  the  variation  of  the  j'th 


rate  coefficient  k,,  clearly  has  a  larger  effect  on  the  concentration  c.(t) 

J  1 


than  the  variation  of  the  other  rate  coefficients.  If  on  the  other  hand  all 


the  coefficients  Av  ',  i  =  l,2,...,n,  are  of  the  same  order  or  magnitude, 
then  the  concentration  of  species  i  at  time  t,  c.(t),  is  effected  essentially 
equally  by  the  variation  of  any  of  the  rate  coefficients  L . 

One  probl em  with  the  above  analysis  must  be  pointed  out.  The  Fourier 
coefficient  of  Eq.  (2.18)  may  be  small  either  because  8c. /Du  is  small  over 
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the  whole  range  of  integration  of  Eq.  (2.19),  which  is  the  case  discussed  above, 
or  because  Dc^/Du^.  changes  sign  one  or  more  times  in  the  range  of  integration. 

Thus,  a  small  value  of  ^  does  not  necessarily  imply  that  the  concentration  c. 

&  1 
is  insensitive  to  changes  in  the  rate  coefficients.  The  remedy  for  this  diffi¬ 
culty  would  be  to  find  a  relation  between  the  Fourier  coefficients,  or  a  combination 
of  Fourier  coefficients,  and  an  everywhere  positive  (or  everywhere  negative) 
function  of  the  rate  of  ciiange  of  ttie  concentrations  with  rate  coefficients,  such 
as  for  instance  ((dc./du^)  >.  We  have,  however,  not  been  able  to  establish 
such  a  relationship  and  it  seems  doubtful  that  a  simple  relationship  of  this  form 
exists.  In  carrying  out  computer  calculations  it  should,  however,  not  be  too 
difficult  to  verify  whether  Sc^./Du^  at  any  given  time  t  is  monotone  or  not  in 
the  range  of  integration  over  the  u  space. 

Thus  while  one  can  assert  that  a  large  implies  high  sensitivity  of  the 

j l 

concentration  of  species  i  with  respect  to  chanaes  in  the  rate  coefficient  k„. 
the  converse  statement  does  not  necessarily  follow  without  checking  on  the  mono- 
tonicity  of  dCj/du  as  discussed  above. 
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III.  Choice  of  the  Frequencies  w... 


As  we  have  mentioned  aoove,  the  application  of  Weyl's  theorem  in  going  from  s 
space,  Eq.  (2.6),  to  the  u  space,  Eq.  (2.16),  via  the  (5  space,  Eq.  (2.8),  must  lead  to 
an  error  in  the  analysis  since  we  use  commensurate  (integer)  frequencies.  This 
error  can  be  minimized  by  a  judicious  choice  of  the  integer  frequencies  u.  and 
the  number  and  magnitude  of  the  m  values  of  the  parameter  s  given  by  Eq.  (2.24). 

We  limit  ourselves  here  to  some  qualitative  remarks  which  will,  however,  clearly 
indicate  the  nature  of  tiie  problem. 

The  integer  frequencies  u.  lead,  according  to  Eqs.  (2.7)  and  (2.24),  to  a 
covering  of  the  n  dimensional  §  space  by  an  array  of  points  as  q  takes  on  its 
integer  values  q=l,2,..,m.  Clearly  one  obtains  a  better  coverage  of  the  6  space 
and  thus  reduces  the  error  in  applying  Weyl's  theorem  if  one  can  increase  the 
density  of  points  in  the  n-dimensional  hypercube  and  if  one  can  distribute  the 
points  uniformly  within  the  space.  Since  Z  and  q  are  both  integer,  the  points 
2.s  (mod  2 it )  will  form  a  regular  lattice  in  9  space.  Our  objective  will  be  to 
make  this  point  lattice,  which  is  completely  generated  by  a  unit  cell,  as  uniform 
as  possiole  in  all  n  directions  by  a  judicious  ctioice  of  the  u's.  As  our  measure 
of  uniformity  we  take  a  nypercubic  unit  cell .  Without  further  information  about 
the  behavior  of  the  output  function  this  seems  the  most  reasonable  choice. 

For  a  fixed  number  m  of  s  points,  chosen  for  computational  convenience, 
we  will  try  to  find  a  vector^  (of  the  infinite  set  of  u's)  which  gives  rise 
to  a  nypercubic  unit  cell.  Once  having  done  this  we  can  assert,  without  loss 

nv:y 

of  generality,  that  ~ -  lies  along  one  edge  of  this  'nypercubic  unit  cell  wtiose 
length  is  — -—-L  •  There  are  now  two  ways  to  compute  tne  volume  V  of  toe 


unit  cell.  In  the  n  dimensional  '6  space,  it  follows  from  the  above  construction 


r>l  v  n 
10)1  \ 


(3.1). 


But.  we  also  know  that  the  total  volume  of  the  n-dimensional  6  space  is  (2")n. 
Since  there  are  m  unit  cells  in  that  space,  it  follows  that 


V=I(2»)n 


Equating  (3.1)  and  (3.2)  then  yields 


(3.2) 


(3.3) 


for  the  'relation  between  the  minimum  length  of  tne  vector  w,  the  dimension  n 
and  tne  number  of  s  points,  m,  for  a  point  lattice  composed  of  hypercibic 
unit  cells. 

The  condition  expressed  in  equation  (3.3)  yields  important  insight  into 
the  judicious  choice  of  the  frequencies  u^.  For  systems  with  large  dimensions, 
i.e.  a  large  number  n  of  rate  coefficients  (or  coupling  parameters  in  general), 

M  approaches  m.  Thus,  the  choice  of  the  number  m  of  s  points  for  the  numerical 
computation  determines  the  value  of  u  and  thus  guides  one  in  the  choice  of 
the  u.'s.  Since  one  would  expect  the  error  in  the  analysis  due  to  the  use  of 
integer  frequencies  w.  to  be  of  order  1/m,  i.e.  inversely  proportional  to  the  number 
of  unit  cells,  it  is  evident  that  one  should  cnoose  a  large  value  of  m  to  carry 
out  the  calculations.  This  in  turn  implies  from  Eq.  (3.3)  a  large  |£|. 

The  analysis  leading  to  Eq.  (3.3)  is  based  on  the  construction  of  unit 
cells  which  are  exact  hypercubes.  Since  |J| ,  m  and  n  are  all  integers  it  may 
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not  be  possible  to  fulfill  condition  (3.3)  exactly  for  all  (arbitrary)  choices 
of  id  and  n.  It  is  easy  to  show  this.  Let  us,  for  instance,  square  both  sides 
of  Eq.  (3.3)  and  let  n  =  5.  We  are  then  required  to  find  an  integer  m  such  that 
we  can  express  a  sum  of  n  squares,  |w|2,  as  ni^0.  Clearly,  this  is  possible, 
if  at  al  lt  only  for  certain  special  values  of  m.  To  be  more  realistic,  we  should 
weaken  our  criteria  for  uniformity  of  the  point  lattice  by  stipulating  unit 
cells  which  are  as  close  as  possible  to  hypercubic  and  then  rewrite  Eq.  (3.3)  as 

Hll  (3.4) 

|u|  *  m  n 

For  n  »  1,  for  which  both  (3.3)  and  (3.4)  reduce  to  |£|  =  m,  it  should  be 
possible  to  obtain  a  more  nearly  exact  hypercubic  unit  cell. 

The  question  as  to  the  "optimum"  choice  of  the  u's  has  been  considered  by 
a  number  of  authors,  using  a  different  approacn  from  tnat  presented  above,  in 
connection  with  the  general  problem  of  the  approximate  evaluation  of  multi¬ 
dimensional  integrals  via  discrete  summations^.  Korobov's  book  has  tables 
of  w^'s  for  a  given  number  of  points  m  (for  m  prime)  and  dimension  n  up  to 
n  =  10.  These  tables  are  reprinted  in  the  book  by  A.  H.  Stroud.  It  is  inter¬ 
esting  and  comforting  to  note  that  although  these  tables  were  computed  from 
completely  different  criteria  than  those  employed  by  us,  the  Korobov  u's  indeed 
generate  hypercubic  unit  cells  to  a  very  good  approximation  for  all  his  sets  oi 
for  which  we  have  carried  out  the  appropriate  calculations. 

Korobov's  analysis  of  the  error  in  the  use  of  integer  frequencies  yields 
explicit  prescriptions  for  calculating  the  "optimum"  set  of  w.'s  as  given  in  his 
tables.  Our  approacn  presented  above  does  not  yield  such  an  explicit  algorithm. 
We  plan  to  develop  sucii  an  algorithm  and  then  compare  our  sets  of  Z  with  those 
of  Korobov  in  a  subsequent  publication. 
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frequency 

As  stated  above,  if  one  wishes  to  decrease  the  error  in  the  integer  A  analys 
one  should  use  a  large  number  m  of  s  points  and  thus  a  large  [w|.  It  can 


readily  be  verified  from  Eq.  (2.6)  that  for  large  one  needs  to  evaluate 

c ( s )  for  a  larger  number  of  s  values  in  order  to  obtain  an  accurate  value  for 
the  Fourier  coefficient  A^.  This  requires  more  extensive  computer  calculations 
A  reasonable  compromise  between  these  two  effects  needs  to  be  adopted 

The  transformation  from  u  space  to  3  space  as  given  by  Eqs.  (2.3)  and  (2.7) 
will  also  effect  the  error  term  since  the  specific  transformation  which  is  chosen 
determines  the  rate  of  change  of  the  function  c(s)  as  a  function  of  s.  We  are 
faced  nere  with  an  interesting  problem  in  compensating  effects.  Either  the 
transformation  u^  =  u^ (Q)  is  singular  or  the  weight  function  p(u^ ; a . )  is  singu¬ 
lar  at  one  or  more  values  of  u^ .  For  the  chosen  transform  it  is  readily  veri¬ 
fied  from  Eqs.  (2.15)  and  (2.20)  that  u.  diverges  at  9  =  n/2  and  6  =  3;i/2,  but 
at  these  points  p(u.;a^)  =  0.  This  same  effect  will  be  found  for  any  transfor¬ 
mation  and  its  associated  weight  function.  Thus  in  the  regions  of  0  space 
where  the  transformation  is  divergent  the  associated  weight  function  will 
always  compensate.  We  are  therefore  led  to  expect  that  the  choice  of  the 
transformation  function  will  not  significantly  effect  the  final  numerical  results. 
This  is  born  out  by  the  data  presented  in  paper  II. 
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ABSTRACT 


The  Fourier  amplitude  method  developed  in  Part  I  as  a  diagnostic  tool 
for  determining  the  sensitivity  of  the  results  of  complex  calculations  to  the 
parameters  which  enter  these  calculations  has  been  applied  to  two  chemical 
reaction  systems  involving  sets  of  coupled,  non  linear  rate  equations.  These 
were:  a)  a  5  reaction  set  describing  the  high  temperature  (6000°K)  dissociation 
of  aii  and  b)  a  9  reaction  set  describing  the  constant  temperature  (2000°K)  com¬ 
bustion  of  H2  and  Og.  We  have  evaluated  the  Fourier  amplitudes  for  all  the 
species  at  a  number  of  different  times  for  both  reaction  systems.  The  analysis 
of  these  results  verifies  the  claims  made  in  part  I.  The  relative  magnitudes 
of  the  Fourier  amplitudes  showed  a  several  order  of  magnitude  distribution  which 
permitted  a  clear  distinction  of  the  relative  sensitivity  of  the  species  concen¬ 
tration  to  uncertainties  in  the  rate  coefficients.  The  conclusions  based  on  the 
Fourier  amplitude  method  for  these  two  reaction  systems  are  in  excellent  agreement 
with  sensitivity  predictions  which  could  be  made  on  the  basis  of  previous  kinetic 
studies  of  these  systems. 


I.  INTRODUCTION 


In  the  preceding  paper,  hereafter  referred  to  as  we  developed  a 
diagnostic  technique  to  investigate  the  sensitivity  of  the  solution  of  large 
sets  of  coupled,  non-linear  equations  to  uncertainties  or  changes  in  the  coup! 
parameters.  For  brevity,  we  will  refer  to  this  method  as  FAST,  the  Fourier 
Amplitude  Sensitivity  Test.  In  this  paper  we  wish  to  demonstrate  the  utility 
of  FAST  by  applying  it  to  a  study  of  the  sensitivity  of  two  sets  of  coupled 
non-linear  rate  equations  describing,  respectively,  the  high  temperature 
dissociation  of  air  and  the  high  temperature  IU  -  0^  reaction.  The  specific 
objective  of  the  calculations  described  here  was  to  determine  the  sensitivity 
of  the  various  species  concentrations  at  a  particular  time  to  assumed  uncer¬ 
tainties  in  the  rate  coefficients  which  enter  into  these  calculations. 

In  Section  II  we  summarize  the  steps  which  must  be  carried  out  in  the 
application  of  FAST  to  any  system  of  coupled  non-linear  equations  and  discuss 
the  choice  of  frequencies  and  the  number  of  sampling  points  to  be  used.  In 
Section  III,  we  present  the  sensitivity  analysis  for  the  two  chemical  reaction 
systems.  Section  IV  contains  some  concluding  remarks. 

II.  GENERAL  REMARKS  ON  APPLICATION  OF  FAST 

We  wish  to  determine  the  sensitivity  of  an  "output  function",  in  our 
case  the  calculated  concentration  of  a  particular  species  at  a  particular  time 
to  assumed  uncertainties  in  the  rate  coefficients  entering  the  rate  equations. 
We  write  this  output  function  as 

c ( ^1  *  ^2  *  '  ‘  ’  ’  ^n ^  c  ( ^)  (2.1) 

where  t  is  the  n-  component  vector  of  the  rate  coefficients  k-j ,  ...,  kf). 


The  application  of  FAST  is  carried  out  via  the  following  steps; 

Step  1)  One  must  choose  a  set  of  n  integer  frequencies,  {u^,  o^,  u  }, 
where  n  is  the  number  of  independent  rate  coefficients  which  determine  the 
value  of  the  output  function. 

Step  2)  One  frequency  of  this  set  is  now  assigned  arbitrarily  to  each 
rate  coefficient  by  (see  1.2.2): 


k,  -  kj°>e 
u.  =  uJo)  sin  u -s 


(2.2) 


(2.3) 


v.hei  e  k.j  is  the  best,  estimate  of  the  rate  coefficient,  s  is  a  parameter  and 
the  u.(o)  . 

i  are  the  endpoints  of  the  estimated  ranges  of  uncertainties  of  the 
rate  coefficients  .  If  a  rate  coefficient  is  known  precisely,  i.e.  with  zero 
uncertainty,  then  k.  =  k!o)  which  implies  u!o)  -  0.  The  endpoints  u!°>  are 
specified  as  part  of  the  input  data  and  their  values  should  be  based  on  the  best 
judgement  of  the  investigator. 

Step  3)  The  output  function  c(£)  =  c[lc(s)]  is  now  Fourier  analyzed  in 
s  to  obtain  the  sine  amplitudes  A^_  corresponding  to  the  input  frequencies  u . . 
This  analysis  requires  the  evaluation  of  c [1c ( s ) ]  for  N  values  of  s,  where  N  is 
an  integer  which  depends  on  the  frequency  spectrum  of  c[£(s)]. 

These  steps  will  now  be  discussed  in  more  detail, 
a)  Choice  of  Frequencies  uk 

In  section  III  of  part  I  we  have  presented  a  discussion  of  the  choice 
of  the  set  of  frequencies  {w}  and  have  also  given  there  a  set  of  references 
(ref.  4)  which  contain  both  discussions  of  optimal  choices  of  sets  of  frequencies 
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as  well  as  tables  of  frequencies  which  could  be  used  for  FAST.  For  the  calcu¬ 
lations  presented  here  we  have  developed  our  own  set  of  frequencies.  We  have 
done  so  to  obtain  a  better  insight  in  how  to  handle  the  calculations  involved 
in  the  Fourier  analysis  method.  In  a  follow-on  study  we  plan  to  investigate 
whether  our  method  of  choosing  frequencies  yields  better,  equivalent  or  worse 
sets  for  use  in  FAST  than  those  presented  in  reference  4  of  part  I. 

Since  we.  will  carry  out  computer  calculations,  we  are  limited  to  rational 
(or,  equivalently,  integer)  frequencies.  For  the  coupled  non-linear  equations 
under  study  here,  this  means  that  in  addition  to  the  input  frequencies  w.  and 
their  multiples,  various  linear  combinations  of  integer  multiples  of  w. 
("interference  frequencies")  will  appear  in  the  spectrum  of  c[1c(s)].  This 
presents. a  problem  in  our  analysis  since  the  Fourier  sine  amplitude  of  frequency 
“a  is  to  reflect  only  the  effect  of  the  uncertainty  in  the  value  of  the  rate 
coefficient  kp.  Clearly,  if  an  interference  frequency  coincides  with  one  of  the 
input  frequencies  to. ,  say  w; ,  the  corresponding  Fourier  amplitude  A  will  not 
only  reflect  changes  in  k}  but  also  in  other  rate  coefficients.  Thus,  for 
example,  if  w,  +  co„  -  w.  =  w9,  the  amplitude  A  is  identical  with  A,  ,  , 

i  J  h  <-  w 2  (to-l+tog— to^) 

and  will  therefore  reflect  not  only  changes  in  k2,  but  also  in  ,  k3  and  k^. 

One  therefore  needs  to  choose  a  linearly  independent  set  of  input  fre¬ 
quencies  ok  such  that 


for 


L  aiwi  *  0  ; 


1  = 


a.  integer 


(2.4) 


t  +  1 

i-1 


(2.5) 
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v.'here  M  is  a  parameter  at  the  disposition  of  the  investigator.  We  refer  to 
such  a  set  as  being  free  of  interferences  to  order  M.  The  larger  the  chosen 
value  of  M,  the  greater  the  likelihood  that  the  Fourier  amplitude  of  each  input, 
frequency  reflects  solely  the  uncertainty  of  the  corresponding  rate  coefficient. 
On  the  other  hand,  as  discussed  in  section  III  of  part  I  and  below,  the  larger 


the  chosen  value  of  M,  the  larger  the  maximum  value,  w  ,  of  the  input  f re- 

max  1 


quencies  wnich  will  still  satisfy  (2.4)  and  (2.5)  and,  correspondingly,  the 
larger  the  set  N  of  points  s  which  will  be  required  for  the  evaluation  of  the 
Fourier  amplitudes.  For  the  calculations  presented  in  this  paper  we  have  chosen 
M  ~  4.  In  Table  I  we  present  sets  of  frequencies  which  are  free  of  interferences 
to  4th  order  for  systems  with  dimensionality  n  from  5  to  19,  i.e.  for  systems 
with  from  5  to  19  independent  rate  coefficients.  These  frequencies  were 
determined  via  computer  by  a  trial  and  error  procedure.  The  sets  of  frequencies 


listed  here  have  the  smallest  value  of  u>max  satisfying  conditions  (2.4)  and 


,  un}  are  assigned  to  the  rate  coefficients 


(2.5)  and  are  referred  to  as  "minimal  sets", 
b )  Assignment  of  Frequencies  to  Rate  Coefficients 
The  frequencies  , 

{k-|,  k£,  ...,  kn)  arbitrarily  since  the  results  of  these  calculations,  in  order 
to  be  useful,  must  be  independent  of  both  the  frequency  set  and  the  assignment 
of  the  frequencies.  To  check  that  the  spectrum  of  the  output  function  does  not 
contain  interfering  frequencies  which  coincide  with  the  input  frequencies,  and 
thus  affect  the  results,  one  should  reassign  the  frequencies  to  the  various  rate 
coefficients  and  repeat  the  Fourier  analysis.  If  the  calculated  Fourier  ampli¬ 
tudes  are  invariant  in  magnitude  and  sign  to  these  permutations,  then  the  set  of 
frequencies  is  free  of  interferences  to  the  order  of  the  dimensionality  of  the 
set  of  rate  equations  and  the  results  of  these  calculations  are  reliable.  If, 
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on  the  other  hand,  some  of  the  Fourier  amplitudes  vary  greatly  with  the  pennu- 
tions  of  the  frequency  assignments,  a  better  set  of  input  frequencies  must  be 
generated.  By  "better"  we  refer  to  a  set  {WJ  with  a  larger  parameter  M  and  a 
h  1 9 h e r  a’IT1  c x •  WeJ] ave  carried  out  such  a  permutation  anal ysis  for  the  H  -  0 ? 
sy.-.'te.n.d i sctissed  below  and  verified  tlia t  cur  sets  of  input  frequencies  yielded 
EMiqLa'npl jtudgs.  which  were  invariant  to  the  frequency  jicrmutations . 


c)  Calculation  of  Fourier  Anpli tudes 

As  discussed  in  section  II  of  part  I,  we  are  only  interested  in  the 
Fourier  sine  amplitudes.  For  an  input  frequency  w  ,  the  amplitude  is 

•  2tt 


A 


rC.  TT 

J  c[h(s )]  sin  w£s  ds 


(2.G) 


To  compute  this  integral,  one  must  evaluate  the  function  c [t(s)]  at  a  set  of 
Points  in  the  interval  0  s;  s  £  27t.  In  the  absence  of  any  information  on  the 
form  of  the  output  function  c  we  take  the  points  to  be  equally  spaced  in  that 
interval , 


~  ~  »  9  =  ^  >  2 ,  . . . ,  N 


(2.7) 


where  N  is  some  integer.  With  these  equally  spaced  points  we  can  now  conveniently 
relate  the  Fourier  coefficient  A  as  defined  by  the  integral  (2.6)  to  the  sum 


N 


A‘  =  n-  £  sin  s  c[«(sj] 


q=i 


t  q 


(2.8) 


wliere 


A* 

O), 


=  A 


(2.9) 
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The  error  term  c  is  given  by 


rr> 

V"*  y' 

i  fmm* td 


A(m!!-to . ) 


£  111- 1  j 


(2.10) 


where  Uj  is  any  frequency  in  the  spectrum  of  /[l;(s)]  which  satisfies 


mN  -  u-  =  u 

J  K 


(2.11) 


lliis  error  term  is  duo  to  "aliasing"  which  inevitably  occurs  whenever  a  finite 

number  of  points  are  chosen  on  an  interval  to  evaluate  Fourier  amplitudes.  For 

example,  if  N  equally  spaced  points  are  used,  the  amplitudes  for  w  will 

«£  «- 

unavoidably  include  the  amplitudes  of  a  component  of  frequency  w .  present  in 

J 

cis)  which  satisfies  the  relation  =mN-on .  where  m  is  an  arbitrary  integer. 

The  value  of  N  is  chosen  to  be  the  smallest  integer  such  that 


E  Vi  1  mN 

i-1 


b .  ,m  integers 


(2.12) 


E  ibi’  "H  +  1 

i-1 


(2.13) 


The  Fourier  amplitudes  A  have  been  evaluated  using  Fqs.  (2.8)  through  (2.10). 

£ 

If  the  output  function  c(k)  could  be  expressed  as  a  polynomial  of  degree 
less  than  or  equal  to  M  in  the  variables  u. ,  i  =  1,  ...,  n,  where  (see  Eq.  2.2) 


i  (o) 

Uj  •  In  (H./k  ) 

J  3  g 

then  the  conditions  (2.12)  and  (2.13)  imply  that  r  =  0  so  that 

'W£ 

A*  -  A 

10,  c.)  .■  • 

t  £ 


(2. 1  A) 


(2.15) 
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TIig  conditions  (2.11)  through  (2.13)  also  cleteruinc  \\  the  minimum  number  of  s 
points  required  for  the  evaluation  of  c [l\ ( s ) ] .  These  are  listed  in  fable  1. 

The  Fourier  amplitudes  are  related  to  a  multiple  integral  of  an  output 
function  over  the  uncertainty  space  (see  I.  2.12).  One  way  of  obtaining  the 
Fourier  amplitudes  without  the  troublesome  interferences  and  aliasing  involved 
in  the  me  thou  discussed  above  would  be  to  evaluate  ttie  n-diniensional  integral 
of  Fq .  (I.  2.12)  directly.  He  have  done  this  via  simple  Monte  Carlo  integrations. 
Our  results  show  that  the  number  of  integration  points  required  by  the  Monte  Carlo 
method  to  give  a  comparable  accuracy  is  much  greater  than  the  number  required 
in  the  Fourier  analysis.  Stated  another  way,  for  a  given  number  of  integration 
points  (we  used  200),  the  degree  of  accuracy  of  the  Fourier  method  detailed 
above  far  surpasses  that  of  the  Monte  Carlo  integration, 
d )  Choice  of  the  l.'eig ht  Func tion 

As  emphasized  in  I,  our  final  result,  i.e.  the  relation 

07  <Dci/Du£>  (2.1G) 

involves  a  weight  function  p(u.)  which  can  be  interpreted  as  a  distribution 
function  in  u  space  which  weights  the  uncertainty  in  the  rate  coefficients.  In 
our  treatment  in  part  I,  we  used  the  weight  function  (see  Fq.  2.20) 

a . 

^ ^ 11  i  ^  cosh  a  . u 7  (2.17) 

which  had  several  desirable  properties.  It  is,  for  instance,  bell  shaped  about 

u --0  corresponding  to  k.-kj°\  the  "best"  value  of  the  rate  coefficient  k.. 

i  it  i 

For  computational  convenience,  we  have  used  two  different  weight  functions  for 
the  calculations  presented  in  section  III  below  corresponding  to  two  different 
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transformations  u.  =  u(w^s).  As.  we  had  hoped  and  expected,  our  results  did  not 
depend  significantly  on  the  choice  of  the  weight  function.  This  indicates  that 
our  output  functions  c(k)  were  "well  behaved",  i.c.  did  not  have  any  large  local 
excursions  in  u  space.  It  is  also  in  agreement  with  the  discussion  in  the  last 
paragraph  of  section  III  of  part  I. 

As  is  shown  in  port  I,  the  form  of  the  weight  function  p(u.)  depends  upon 
the  choice  of  the  transformation  function  f  of  Fq.  (1.2.3).  For  the  Fourier 
amplitude  method  to  be  at  all  useful,  our  results  on  the  relative  magnitudes  of 
the  Fourier  amplitudes  must  clearly  be  independent  of  the  explicit  form  of  the 
transformation  function.  The  insensitivity  of  our  results  to  our  choice  of 
transformation  function  (and  thus  of  the  weight  functions)  verifies  the  validity 
of  FAST  in  this  respect. 

III.  APPLICATION  OF  FAST  TO  CHEMICAL  SYSTEMS 

We  have  applied  FAST  to  two  relatively  simple  reaction  systems  as  a 
test  of  its  utility  as  a  diagnostic  tool.  Both  systems  were  assumed,  for 
simplicity,  to  react  at  constant  temperatures  and  volume.  These  restrictions 
are  in  no  way  necessary  for  the  application  of  FAST.  The  concentrations  were 
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integrated  in  time  with  a  code  using  Gear's^  algorithm.  These  integrations 
wet  e  carried  out  once  for  each  of  the  II  value.'  of  s  with  the  corresponding 
values  of  the  n  rate  coefficients,  'fc(s),  up  to  some  time  t.  This  yielded  the 

output  function  c^[k(s),t],  i.e.  the  concentration  of  species  i,  i~l,  2 . m, 

as  a  function  of  time.  The  Fourier  analysis  of  c. [!"($), t]  was  then  carried  out 
at  several  selected  times.  The  uncertainties  in  the  rate  coefficients  were 
arbitrarily  assumed  to  be  much  larger  than  warranted  by  the  available  experi¬ 
mental  data  in  order  to  check  out  FAST.  It  should  therefore  be  pointed  out 
again  that  the  primary  purpose  of  the  calculations  presented  here  is  to  demonstrate 
the  validity  and  utility  of  FAST  as  a  diagnostic  tool  and  not  to  provide  new 
information  about  t fie  selected  reaction  systems. 
a )  High  Temperature  Air  React io n 

He  consider  the  following  simplified  reaction  system'3^ 


09  =  1-1  T— < 

u 

20  +  M 

kl  *  K1 

(1) 

N2  h-  M  — 

2N  +  H 

k2  >  l->2 

(2) 

NO  +  M 

N  +  0  + 

M  k3,  K3 

(3) 

N2  +  0 

NO  +  N 

k4  >  ^4 

(4) 

02  +  N 

NO  +  0 

k5’  K5 

(5) 

The  rate  coefficients  k. 

i 

,  1=1, 

. ,  5 ,  refer  to 

the  forward 

reactions,  the  rate 

coefficients  k  ^  for  the 

reverse 

reactions  are 

obtained  from 

the  equilibrium 

constant  K.  through'  the  relation  Ih  =  k./k_..  In  the  system  and  the  one  dis- 

cussc  below  in  1 1 1 b ,  the  equilibrium  constants  are  assumed  to  be  known  exactly, 
i.e.  with  zero  uncertainty.  The  forward  and  reverse  rate  coefficients  k . ,  k  . 
for  each  reaction  therefore  have  identical  uncertainties  and  the  number  of 
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i ml c pendent  uncertainties  is  equal  to  the  number  of  equilibrium  constants.  It 
is,  of  course,  not  necessary  to  proceed  this  way  and  one  can  instead  assign 
di Ffci  ent  uncertainties  to  the  i onward  and  reverse  rate  coefficients  which  reflect 
uncertainties  in  the  equilibrium  constants. 

The  system  was  assumed  to  react  at  a  constant  temperature  of  6000°K  and 
constant  volume.  The  rate  coefficients,  equilibrium  constants  and  initial  con¬ 
ditions  used  in  this  calculation  are  shown  in  Table  II. 


The  equilibrium  constants  represent  the  latest  and  "best"  literature  values  we 
could  find  and  the  rate  coefficients  listed  in  this  table  are  the  k!0^  of 
Eq .  (2.2).  In  Mg.  1  we  show  a  plot  of  the  time  evolution  of  this  reaction 
system  obtained  from  the  integration  of  the  rate  equations  using  the  rate 
coefficients  k '  of  Table  II. 

While  the  real  uncertainty  of  these  rate  coefficients  is  riot  too  large, 
we  assumed  arbitrarily  for  the  purpose  of  testing  FAST  that  the  experimental 


uncertainties  were  ~  two  orders  of  magnitude  for  each  of  the  five  independent 
sets  of  rate  coefficients.  To  take  account  of  this  uncertainty  we  then  write 
(see  Eq.  (2.2)) 

k(°)  e-4.606  ^  ^  k(o)  e4.606 


i  =  1 . 5 


(3.1) 


For  our  sensitivity  analysis  we  used  the  frequency  set  (see  Tabic  I) 


{to}  =  (2,  42,  62,  74,  90) 


(3.2) 


for  which  the  smallest  number  N  of  evenly  spaced  points  satisfying  the  conditions 
of  Eqs .  (2.4,  2.6,  2.12,  2.13)  in  the  s  interval,  Ores  re  2n,  is  191. 

We  carried  out  the  Fourier  analysis  of  the  output  function  c.[lt(s)]  at 


t  =  10  ^  sec  and  lO-1^  sec.. 


From  Eqs.  (2.2),  (2.3)  and  (3.1)  we  have 


u.j  =  4.606  sin  w- s 
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The  concentrations  of  the  various  species,  ! 12>  0£ >  NO,  0,  I!,  were  calculated 
at  10"6  sec  and  10“4  sec  for  each  of  the  191  s  values  by  integrating  the  rate 
equations .  The  191  values  of  the  concentration  of  the  above  species  at  a  given 
time  define  the  output  function  c.jl.'k(s)]  which  is  now  Fourier  analysed  by 
computing  the  Fourier  coefficients  from  Fqs.  (2.8)  through  (2. Id)  using  the  set 
of  input  frequencies  of  Eq.  (3.2).  An  example  of  the  function  c^|.T.(s)]  for  MO 
at  t  -  10"4  sec  is  shown  in  Fig.  2.  This  curve  has  been  drawn  by  connecting 
the  191  calculated  concentrations  of  MO  (one  for  each  of  the  191  evenly  spaced 
values  of  s)  by  straight  lines. 

The  results  of  the  Tourier  analysis  at  t  -  10-6  sec  are  shown  in  Table 
III.  The  high  values  of  the  Fourier  amplitudes  associated  with  the  rate  coef¬ 
ficients  k-j ,  k_.j  for  the  concentration  of  0  and  0^  (see  blocks  1  and  2  of 
Table  III)  indicate  that  a  variation  in  k-j  and  k_1  lias  a  greater  effect  on  the 
concentration  of  0  and  0^  than  the  variation  of  any  of  the  other  rate  coef 1 icioncs 
This  follows  immediately  from  Eq.  (2.16),  where  now  i  refers  to  the  0  atoms 
and  02  molecules  and  =  2  refers  to  the  rate  coefficients  k-j  and  k_-j .  Me  can 
rephrase  this  by  saying  that,  at  this  early  time  in  the  reaction,  only  reaction 
(1)  is  important  in  the  production  of  0  atoms  and  the  disappearance  of  02 
molecules.  The  sign  of  the  amplitudes  in  the  last  column  can  be  understood  from 
Eq.  (2.16)  according  to  which  the  amplitude  A^  is  proportional  to  (dc./Du&>. 
Thus  the  negative  sign  of  the  Fourier  amplitude  for  k-j  for  the  concentration  of 
02  (block  1,  Table  III)  simply  indicates  that  [02]  decreases  with  increasing  k] . 
The  positive  sign  for  the  Fourier  amplitude  for  k-j  for  the  concentration  of  0 
atoms  (block  2,  Tabic  III)  indicates  that  [0]  increases  with  increasing  k-, . 

From  the  relative  magnitudes  of  the  Fourier  amplitudes  in  blocks  4  and  5  of 
Table  III,  it  will  be  note  I  that  the  formation  of  N  and  MO  at  these  early  times 
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in  the  reaction  is  most  sensitive  to  reactions  (1)  and  (4).  This  is  quite 
reasonable  since  at  t  ~  10  ^  sec  the  formation  of  II  and  HO  via  reactions  (?) 
and  (5)  cannot  play  any  important  role  owing  to  the  small  value  of  k2°^  as 
compared  to  kj0^  (see  Table  II).  A  number  of  other  conclusions  as  to  the 
sensitivity  of  the  concentrations  of  the  various  species  with  resect  to  vari¬ 
ation  of  the  rate  coefficients  can  be  drawn  from  a  study  of  the  results  in  Table 
III.  All  these  conclusions  are  in  accord  with  one's  knowledge  and  intuition  of 
this  rather  simple  reaction  system. 

In  Table  IV  we  present  our  results  for  t  =  10"^  sec.  At  this  time,  as 
can  be  seen  from  Fig.  1,  the  concentrations  of  reactants  and  products  are, 
except  for  significantly  different  from  their  initial  values.  The  coupling 

between  the  various  reactions  has  come  into  play  much  more  strongly  at  this  later 
time  and  one  would  expect  that  the  influence  of  the  uncertainties  in  the  rate 
coefficients  on  the  concentrations  of  the  various  species  is  more  complex.  This 
can  readily  be  verified  from  the  data  in  Table  IV,  where  there  is  now  no  longer 
such  a  pronounced  order  of  magnitude  effect  in  the  Fourier  amplitudes.  Let  us 
arbitrarily  use  a  factor  of  about  10  to  distinguish  between  the  "relatively 
important"  amplitudes  (i.e.  the  concentration  of  the  species  is  very  sensitive 
to  the  change  of  a  specific  rate  coefficient)  and  the  "relatively  unimportant" 
amplitudes  (i.e.  the  concentration  of  the  species  is  net.  very  sensitive  to  the 
change  of  a  specific  rate  coefficient).  Then  the  data  of  Table  IV  indicates 
that  [02]  and  [0]  are  affected  primarily  by  uncertainties  in  1^  ,  [N^l  is 
affected  about  equally  by  uncertainties  in  all  of  the  rate  coefficients,  [N]  is 
ir.ost  affected  by  uncertainty  in  k-j ,  and  [HO]  is  least  affected  by  uncertainty 
in  k?.  The  least  "sensitive"  rate  coefficient  is  k2>  which  determines  the  rote 
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of  dissociation  of  N2<  This  is  not  surprising  since  k^°)  is  smaller  than  the 
I  other  rate  coefficients  by  2  to  4  orders  of  magnitude.  The  important  point  to 

note  is  that  owing  to  the  strong  coupling  of  all  the  reactions  at  this  stage  of 
the  kinetic  development,  a  number  of  the  elementary  rate  processes  contribute 
}  significantly  to  the  production  and  removal  of  most  of  the  species,  so  that  the 

concentration  of  these  species  are  quite  sensitive  to  uncertainties  in  several 
rather  than  just  one  rate  coefficient.  This  result  is  certainly  not  unexpected. 

f  b)  The  -  0^  Reaction 

The  H2  -  02  combustion  system  was  modeled  by  the  following  set  of  chemical 
reactions  ^ : 


H2 

+ 

M  r- 

- ^ 

2'!  1 

-  il 

kl  ’ 

h 

(1) 

°2 

+ 

!•!  P-- 

> 

20  i 

■  M 

k2* 

h 

(2) 

OH 

+ 

M  , 

_  > 

0  + 

II  M 

kg, 

K3 

(3) 

h2o 

+ 

H  f= 

II  + 

OH  +  M 

k^ , 

h 

(4) 

H 

+ 

o2<p- 

OH  a 

■  0 

k5‘ 

h 

(5) 

0 

+ 

ll2  i= 

- > 

OH  h 

■  II 

kc> 

h 

(6) 

H2 

+ 

0H  v- 

r-> 

ll20 

+  II 

k7* 

K7 

(7) 

°2 

+ 

h2  f= 

rzf 

2011 

k8, 

K8 

(8) 

OH 

+ 

OH  P 

— } 

H20 

+  0 

> 

K9 

(9) 

The  system  was  assumed  to  react  at  a  constant  temperature  of  2000°K  and  constant 
volume,  with  only  n2  and  02  present  at  t  =  0. 

The  initial  conditions,  rate  coefficients  and  equilibrium  constants  are 
shown  in  Table  V.  The  rate  coefficients  listed  in  this  table  are  the  kH  of 
fq.  (2.2):  the  equilibrium  constants  are  assumed  to  be  known  with  zero  uncertainty 
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We  assumed  arbitrarily  that  the  experimental  uncertainties  for  each  of 
the  nine  independent  rate  coefficients  were  ±  one  order  of  magnitude  so  that 


k(°)  -2.303  .(0)  2. 505 

1  “  *'i  -*  lvi  e 


Fot  our  analysis  we  used  the  frequency  set 


i  -  1 ,  2,  . . . ,  9  . 


(3.4) 


{“)  =  (]9>  59 »  91  »  113,  I33,  143,  149,  1  57,  161  ) 


(3.5) 


for  which  the  smallest  number  II  of  evenly  spaced  points  in  the  s  interval, 

OSS  s2r,  is  630. 

A  plot  of  the  concentrations  of  the  various  species  for  the  time  interval 
10  ^  t  *  10  seconds,  calculated  with  the  l:.[o)  of  Table  V  is  shown  in  Tig.  3. 

We  have  employed  FAST  to  study  the  sensitivity  of  the  species  concentrations  to 
uncertainties  in  the  rate  coefficients  at  10~8  seconds  and  10‘3  seconds.  These 
results  are  shown  in  Tables  VJ  and  VII. 

At  the  very  early  time  of  10"8  seconds  (Table  VI),  the  Hg  and  02  concen¬ 
tration,  have  not  changed  significantly  from  their  original  value  (to  within 
8  significant  places)  and  the  Fourier  amplitudes  are  merely  round-off  error. 

For  H  atoms,  reaction  (1)  seems  to  be  the  major  contributor  to  its  formation 
at  this  early  time  and  its  rate  coefficient  is  the  most  sensitive  in  determining 
the  accuracy  of  [II].  Reactions  (7)  and  (8),  which  together  also  produce  II  atoms, 
are  also  important  but  the  concentration  of  H  atoms  is  not  as  sensitive  to  k? 
and  k8  as  it  is  to  ky  The  sensitivity  of  [H]  to  these  two  rate  coefficients  has 
identical  values  which  is  not  surprising  since  [H]  is  proportional  to  k?-kQ  at 
this  early  time.  The  same  argument  also  holds  for  the  production  of  HO  since 
its  concentration  at  10  Seconds  is  proportional  to  k? •  kQ.  The  ll.o  concentration 
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is  also  sensitive  to  k^,  since  reaction  (Li)  appears  to  be  partly  responsible 
for  the  production  of  the  Oil  which  enters  into  reaction  (7).  The  concentration 
of  0  atoms  is  most  sensitive  to  kg  which  is  again  as  expected.  Finally,  the  Oil 
concentration  is  most  sensitive  to  kg  which  is  also  most  reasonable,  since  at 

this  early  time,  only  the  direct  reaction  between  Og  and  1 1^  would  be  expected 

to  effect  the  Oil  concentration.  Again,  the  results  of  FAST  are  in  good  agreement 
with  what  could  be  predicted  on  the  basis  of  chemical  knowledge, 

_3 

At  t  ••  10  '  seconds  (Table  VII)  the  story  is  quite  different.  As  can  be 

seen  From  Fig.  3,  significant  changes  from  the  initial  concentrations  have 

o ecu red  by  this  time,  and  the  reactions  are  now  strongly  coupled.  This  is 

reflected  in  the  array  of  Fourier  amplitude  in  Table  VII  as  compared  to  those 

of  Table -VI.  How,  for  instance,  rate  coefficient  kn  corresponding  to  the  chain 

branching  step  (5)  which  produces  H  and  0  atoms  is  the  most  influential  one  in 

controlling  the  concentration  of  all  species.  The  II  atoms  concentration,  in 

addition,  is  equally  sensitive  to  k-,.  All  other  rate  coefficients,  except  for 

kg,  which  is  very  small  for  all  species,  have  about  equal  sensitivities  to 

within  a  factor  of  10  as  measured  by  their  Fourier  amplitudes.  According  to 

our  theory  and  calculations  then,  the  rate  coefficient  which  needs  to  be 

determined  with  the  highest  accuracy  for  the  accurate  calculation  of  all  species 

-3 

concentration  at  10  seconds  is  kg  with  k-,  next  in  line.  Uncertainties  in  the 
other  rate  coefficients  will  have  a  smaller  effect  on  the  species  concentrations. 
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IV.  CONCLUDING  REMARKS 


It  should  bo  realized  that  the  applications  presented  here  arc  only  a 
firs'  cut  at  the  subject.  It  appears  that  FAST  is.  a  useful  diagnostic  tool  in 
dotei  mining  the  sensitivity  of  the  results  of  complex  calculations  to  the 
parameters  which  enter  the  calculations.  The  Fourier  amplitude  method  clearly 
need.,  further  testing  on  more  complicated  systems  than  the  ones  dealt  with  in 

this  l>a pc i .  Such  tests  should  also  involve  comparison  with  various  types  of 
"brute  force"  calculations. 

In  addition  to  ap;  Mentions  to  larger  systems,  there  are  also  a  number  of 
theoretical -computational  questions  which  need  further  investigation.  Some 
crumple-  ate.  the  influence  of  the  range  of  uncertainties,  with  different 
uncertainties  for  different  coupling  parameters,  on  the  relative  magnitudes  of 
the  Fourier  amplitudes;  the  influence  of  the  choice  of  the  input  frequencies  (,.,} 
and  the  spacing  of  the  I!  values  of  the  parameter  s  on  the  Fourier  amplitudes;  am 
the  question  whether  the  absolute  magnitude  of  the  Fourier  amplitudes  can  be 
used  as  a  predictor  for  the  actual  change  of  the  output  function  for  a  given 

change  of  the  coupling  parameter.  Additional  work  is  being  carried  out  to  answe. 
these  questions. 
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FIGURE  CART  I OilS 

I'iy*  1  iiiie  evolution  of  the  !{,,  -  0 ^  system  based  on  the  binotic  data 

in  Tab! c  II. 

Fig.  2.  The  function  c.['l;(s)j  for  ilO  at  10"4  sec  between  0  and  2n  as 
a  function  of  s.  This  curve  has  been  drawn  by  connecting  the 
191  calculated  values  of  the  concentrations  of  110  (corresponding 
to  the  191  values  of  s)  by  straight  lines. 

Fig.  3.  Time  evolution  of  the  ll?  -  0^  systems  between  t  -  10~4  sec  to 
t  =  10  sec  based  on  the  kinetic  data  in  Table  V. 
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Table  I 


FREQUENCY  SETS  AND  CORRESPONDING  NUMBER  OF 
POINTS  AVOIDING  INTERFERENCES 
THROUGH  4th  ORDER 


Dimension  (n) 

Frequency  Set  {u} 

Minimum  Number  of  Points  (N) 

5 

11,  21,  27,  35,  39 

142 

5* 

2,  42,  62,  74,  90 

191 

6 

■  1  ,  21  ,  31  ,  37,  45,  49 

182 

G* 

2,  42,  62,  74,  90,  98 

231 

7 

17,  39,  59,  69,  75,  83,  87 

334 

8  • 

23,  55,  77,  97,  107,  113, 

121,  125 

486 

9 

19,  59,  91,  113,  133,  143 

149,  157,  161 

630 

10 

25,  63,  103,  135,  157,  177, 
187,  193,  201,  205 

806 

11 

41,  67,  105,  145,  177,  199, 
219,  229,  235,  243,  247 

974 

12 

31  ,  87,  113,  151  ,  191  ,  223, 
245,  265,  275,  281,  289,  293 

1158 

13 

23,  85,  141,  167,  205,  245 
277,  299,  319,  329,  335,  343, 

1374 

347 

14 

87,  133,  195,  251,  277,  315, 
355,  387,  4U9,  429,  439,  445, 

/l«  /it;7 

1814 

21  -  ^ 


Table  II 


INITIAL  CONDITIONS,  RATE  COEFFICIENTS  AND  EQUILIBRIUM 
CONSTANTS  FOR  THE  HIGH  TEMPERATURE  AIR  REACTION 

Initial  Conditions 
[N2]  =  8  x  10“8  moles/cc 
[02]  =  2  x  10”8 


Rate  Coefficients 


k|0^  =  8.5  xlO^8  (moles/cc)"^  sec  ^ 
=  3.0  x  I07 

k^o)  =  8.0  x  109 

k|o)  =  9.0  x  1010 

kg°)  =  8.0  x  TO11 


Equilibrium  Constants 

K-j  =  7.8  x  10"4  (moles/cc) 

K2  =  1 .0  x  10"7 

K3  =  1.3  x  10"5  u 

K4  =  9.0  x  10'3 

K5  =  5.9  x  101 


*3- 


. . . - . . . _ 


I 

1 

“  ft 

% 

B 

ij 

K 

r 

ft 

r 

FOURIER  AMPLITUDES  OF 
FOR  THE  HIGH 

Table  III 

THE  CONCENTRATIONS  AT  10 
TEMPERATURE  AIR  SYSTEM 

"6  SEC 

3  1 

I 

i 

r 

Input  Frequency  , 

Rate  Coefficient 

Fourier  Amplitude 

1 

5. 

[o2] 

g 

x  10  moles/cc 

|  i 

1 

i 

2 

kl 

-4.47 

* 

j 

r 

90 

k5 

-0.0244 

74 

k4 

-0.0227 

j 

42 

k2 

-4.7  x  10"4 

I 

t 

62 

k3 

-3.9  x  10"4 

* 

»i 

$ 

[0]  x 

109  moles/cc 

I 

t 

2 

kl 

8.53 

■! 

74 

k4 

-0.36 

i 

J 

90 

k5 

0.018 

j 

t 

62 

k3 

-0.029 

42 

k2 

0.0027 

[N2J  X 

g 

10  moles/cc 

f 

1 

2 

kl 

.  -0.389 

j 

•a 

if 

74 

k4 

-0.386 

I 

J 

Table  IV 
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FOURIER  AMPLITUDES  OF  THE  CONCENTRATIONS  AT  10"4  SEC 
FOR  THE  HIGH  TEMPERATURE  AIR  SYSTEM 


Input  Frequency 

Rate  Coefficient 

Fourier  Amplitude 

[O2]  x  10^  moles/cc 

2 

kl 

-1.1 

90 

k5 

-0.147 

74 

k4 

-0.087 

62 

k3 

-0.069 

42 

k2 

-0.044 

[0]  x  10^  moles/cc 

2 

kl 

2.17 

.  62  • 

k3 

0.228 

90 

k5 

0.164 

42 

0.078 

74 

k4 

-0.00146 

g 

[Ng]  x  10  moles/cc 

74 

k4 

-2.6 

2 

kl 

-2.3 

62 

-1.0 

I 
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Table  V 


Initial  Conditions,  the  Rate  Coefficients  and  Equilibrium 
Constants  for  the  H 2  -  0^  Reaction 


Initial  Conditions 


[1^]  =  8.0  x  10"9  moles/cc 
[O2]  =  4.0  x  10”9  moles/cc 


Rate  Coefficients 


-1  -1 


k°  -  5.78  x  104  (moles/cc)”  sec 


k°  =  4.47  x  103 

k°  =  1 .03  x  105 

k°  =  6.79  x  105 

k°  =  1.23  x  1012 

Equilibrium  Constants 


K-|  =  1.59  x  10”  moles/cc 


4.0  x  10^3  (moles/cc)”^sec”^ 


1.51  x  101 


2.20  x  10u 
1 .6  x  lO14 


K2  =  2.69  x  10' 
K3  =  1.14  x  10' 


K4  =  1.63  x  10' 
K5  =  2.37  x  10' 


Rif' 
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Table  VI 


FOURIER  AMPLITUDES  OF  THE  CONCENTRATIONS  AT  10"8  SEC 
FOR  THE  H2  -  02. REACTION 


Rate  Coefficient _ Fourier  Amplitude 

[h2i. 

Insensitive  to  uncertainties 
in  any  rate  constant. 

J9 


<1 


[H]  x  101 


4.67 

0.202 

0.202 

0.010 

0.0023 

0.00165 


8.0x1 0" 


4.1x10 


-5 


-3.0x10" 


[OH]  x  1018 


_Rate  Coefficient _ Fourier  Amp!  itude 


[0,] 


Insensitive  to  uncertainties 
in  any  rate  constant. 

[0]  x  1020 


1.79 
-0.016 
0.0065 
0.0065 
2.25x10 
-1.95x10' 
7.87x10' 
1 .20x10" 
-8.46x10" 


[H?0]  x  10 


,20 


St*7- 


i  '  [•  li  HHMiyiilili 


tMiSuMwnbuJUM 

— ••* . . 


5 

k8 

5.9 

k8 

if 

2.03 

m 

f 

k7 

-0.0203 

k7 

2.03  If 

1 5 

k9 

0.00176 

k5 

0.121  || 

k5 

-.0.0014 

k6 

p  j 

0.0066 

Y 

k2 

1.61x1  O'4 

‘  k4 

A  |1 

8.35x10“^ 

k6 

9.73xl0"5 

kl 

-6. 35x1 0"4 

£  ,, 

I 

kl 

7. 08x1 0"5 

k9 

4.22x1 0"4  || 

S. 

> 

k4 

-7.43xl0"6 

k3 

-2.81xl0"4  | 

k3 

3.13xl0"6 

k2 

1 . 33x1 0“4 

ivsa^ssgij 


-  . 


rnamems^mmin 


Table  VII 


» 

FOURIER  AMPLITUDES  OF  THE 

CONCENTRATIONS  AT 

10"3  SEC 

FOR  THE  H2  ■ 

•  02  REACTION. 

e# 

IT 

Rate  Coefficient 

Fourier  Amplitude 

Rate  Coefficient  Fourier  Amplitude 

CH23 

x  109 

[H]  x  109 

k5 

-3.71 

k5 

1.92 

r 

k7 

-0.346 

k7 

1.92 

k6 

-0.313 

k6 

0.122 

k8 

-0.172 

k8 

0.102 

r  ■ 

k2 

-0.108 

k2 

0.061 

k4 

• 

-0.101 

k9 

0.0466 

k9 

-0.0905 

k4 

0.0439 

r 

kl 

-0.0491 

kl 

0.0176 

k3 

0.0119 

k3 

-0.00318 

Co2] 

X 

o 

VO 

! 

[0]  x  1010 

i 

k5 

-1.77 

k5 

5.77 

k7 

-0.153 

k7 

0.420 

i 

k6 

-0.089 

kl 

0.164 

k8 

-0.070 

k6 

0.162 

9i 
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II.  APPLICATIONS 


1.  R.I.  Cukier,  C.N.  Fortuin,  K.E.  Shuler,  A.G.  Petschek 
and  J.H,  Schaibly ,  J.  Chem.  Phys . (preceding  paper) 


C.W.  Gear,  "The  Numerical  Integration  of  Stiff 
Differential  Equations,"  University  of  Illinois,  Dept, 
of  Computer  Sciences,  Report  No. 221,  Jan.  1967. 


The  equilibrium  constants  were  calculated  from  data  in 
B.J.  McBride,  Thermodynamic  Properties  to  6000°  K 
_for  210  Substances  Involving  the  First  18  Elements. 


NASA,  SP-3001  (1963).  The  rate  coefficients  were 
taken  from:  J.C.  Schexnayder  and  J.S.  Evens,  NASA 
Technical  Report  R-108,  (1961);  H.S.  Click,  J.J.  Klein, 
and  W.  Squire,  J.  Chem.  Phys.,  27_,  850  (1957); 

G.B.  Kistiakowski  and  G.C.  Volpi,  J.  Chem.  Phys. 27, 
1141,  (1957).  We  are  aware  that  more  recent  "nominal" 
values  of  the  various  rate  coefficients  exist  in  the 
literature  but  feel  that  the  values  used  in  our 
calculations  are  adequate  for  the  purpose  of  this 
paper . 


The  equilibrium  constants  were  calculated  from  the 
data  in  McBride  (see  Ref. 2.);  the  rate  coefficients 
were  obtained  from  the  compilation  of  G.S.  Bahn, 
Reaction  Rate  Compilation  for  the  H-O-N  System, 
Gordon  and  Breach  Science  Publishers,  New  York  (1968) 
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